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Abstract. In this article we study the semi-classical distribution of the complex zeros of the eigenfunctions 
of the ID Schrodinger operators for the class of real polynomial potentials of even degree, with fixed energy 
level, E. We show that as /in ^ the zeros tend to concentrate on the union of some level curves 



?R{S{zm, z)) = Cm where S{zm, z) = sJV{€) — E dt is the complex action, and Zm is a turning point. We 
^ ' also calculate these curves for some symmetric and non-symmetric one- well and double- well potentials. The 

example of the non-symmetric double- well potential shows that we can obtain different pictures of complex 
zeros for different subsequences of hn. 

Keywords: Schrodinger operator, Complex WKB method. Stokes lines. 

1. Introduction 

^ ' This article is concerned with the eigenvalue problem for a one-dimensional semi-classical Schrodinger 

. operator 

o 

(1) {-h^-^ + V{x))tlj{x,h) E{h)tl){x,h), ij{x,h)eL\R) h^0+ 

' Using the spectral theory of the Shrodinger operators [BSj, we know that if lim^^^oo — +oo then the 

, spectrum is discrete and can be arranged in an increasing sequence £^o(^) < Ei(h) < E2{h) <■■■] oo. Notice 

' that each eigenvalue has multiplicity one. We let {V'n(a;, /i)} be a sequence of eigenfunctions associated to 

En{h). If we assume the potential V{x) is a real polynomial of even degree with positive leading coefficient, 
' then we can arrange the eigenvalues as above and the eigenfunctions ipnix, h) possess analytic continuations 

. 'i'niz, h) to C. Our interest is in the distribution of complex zeros of ipniz, h) as h —>■ 0+ when an energy 

level E is fixed. The substitutions A = -^j and q{x) = V{x) — E changes the eigenvalue problem ([J) to the 
problem: 

(2) y"{x,\) = \^q{x)y{x,\), A) G ^^(R), A ^ oo. 
Since limj-^+oo 9(2;) = +00, again the spectrum is discrete and can be arranged as 

(3) Ao < Ai < A2 < • • • < A„ • • • T 00. 
We define the discrete measure Z\^^ by 

(4) ^^"^^ E ^- 

" {z\ y(2,A„)=0} 

In this paper we study the limits of weak* convergent subsequences of the sequence {Z\^} as n ^ 00. We 
say 

Z\^^ — > Z, (in weak* sense) 

if for every test function ip G (R^ ) we have 

^A„J'P) — > Z{lp). 

We will call these weak limits, the zero limit measures. 

Throughout this article we assume that q{z) has simple zeros. We may be able to extend the results in the 
case of multiple turning points using the methods in \ P1\ I02j on the asymptotic expansions around multiple 
turning points. Notice that q{x) has to change its sign on the real axis, because if it is positive everywhere 
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then ((2]) does not have any solution in L^(M). Hence q{z) has at least two simple real zeros. We say q{x) is 
a one-well potential if it has exactly two (simple) real zeros and a double- well potential if it has exactly four 
(simple) real zeros. One of our results is 

Theorem 1.1. Let q{z) be a real polynomial of even degree with positive leading coefficient. Then every 
weak limit Z (zero limit measure) of the sequence {2^a„} is of the form 

(5) Z=-\^)\\djl 

TT 

m^W{t)\dt, 

where j is a union of finitely many smooth connected curves 7„ in the plane. For each 7^ there exists a 
constant Cm, « canonical domain D,n and a turning point z^ on the boundary of such that 7^ is given 
by 

(6) ^{S{Zm,z)) = Cm, Z e Dm, 

where S{zm, z) = ^ q(t) dt and the integral is taken along any path in Dm joining Zm to z.(See section 
2, for the definitions) 

This theorem shows that if Z in ([5]) is the limit of a subsequence {Zx^^}, then the complex zeros of 
{y{z, Xn^)} tend to concentrate on 7 as fc — ^■ 00 and in the limit they cover 7. The factor q{z)\ = 
^yV(z) — E\ indicates that the limit distribution of the zeros on 7 is measured by the Agmon metric. We 
call the curves jm the zeros lines of the limit Z. The next question after seeing Theorem 11.11 is "what are 
all the possible zero limit measures and corresponding zero lines for a given polynomial q{z)7" We answer 
this question for some one- well and double- well potentials. But before stating these results let us mention 
some background and motivation for the problem. 

Form the classical Sturm-Liouville theory we know everything about the real zeros of solutions of ^ . We 
know that on a classical interval (i.e. an interval where q{x) < 0), every real-valued solution y{x,X) of Q 
(not necessarily L^-solution) is oscillatory and becomes highly oscillatory as A ^ 00. In fact the spacing 
between the real zeros on a classical interval is j . On the other hand there is at most one real zero on each 
connected forbidden interval where q{x) > 0. This shows that every limit Z in ([5]) has the union of classical 
intervals in its support. 

It turns out that other than the harmonic oscillator q{z) ~ z^ — cP' where the eigenfunctions do not have 
any non-real zeros, the complex zeros are more complicated. It is easy to see that when q(z) = z^ -\- az^ -\- b, 
the eigenfunctions have infinitely many zeros on the imaginary axis. For q{z) = z'^+a'^, Titchmarsh in [T] 
made a conjecture that all the non-real zeros are on the imaginary axis. This conjecture was proved by 
Hille in [Hlj . In general one can only hope to study the asymptotics of large zeros of y{z, A) rather than 
finding the exact locations of zeros. The asymptotics of zeros of solutions to ([2]) for a fixed A and large 
z, have been extensively studied mainly by E. Hille, R. Nevanlinna, H. Wittich and S. Bank (see [N], |Wj . 
[B]). But it seems the semi-classical limit of complex zeros has not been studied in the literature, at least 
not from the perspective that was mentioned in Theorem ll.il which is closely related to the quantum limits 
of eigenfunctions. This problem was raised around fifteen years ago when physicists were trying to find a 
connection between eigenfunctions of quantum systems and the dynamics of the classical system. It was 
noticed that for the ergodic case the complex zeros tend to distribute uniformly in the phase space but for the 
integrable systems the zeros tend to concentrate on one- dimensional lines. An article which made this point 
and contains very interesting graphics is |LVj . The problem of complex zeros of complexified eigenfunctions 
and relations to quantum limits is suggested by S. Zelditch mainly in [Z]. There, the author proves that if a 
sequence {(p\„ } of eigenfunctions of the Laplace-Beltrami operator on a real analytic manifold M is quantum 
ergodic then the sequence {•Z'a„} of zero distributions associated to the complexified eigenfunctions {tp^ } on 
M^, the complexification of M, is weakly convergent to an explicitly calculable measure. A natural problem 
is to generalize the results in [Z] for Schrodinger eigenfunctions on real analytic manifolds. This is indeed 
a difficult problem. Perhaps the first step to study such a problem is to consider the one dimensional case 
which we do in this paper. The main reason to study the complex zeros rather than the real zeros is that the 
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problem is much easier in this case (in higher dimensions) . For example in studying the zeros of polynomials 
as a model for eigenfunctions, the Fundamental Theorem of Algebra and Hilbert's NuUstellensatz are two 
good examples of how the complex zeros are easier and somehow richer. See [Zlj for some background of 
the problem and some motivation in higher dimensions. 

One of our results is that for a symmetric quartic oscillator the full sequence {Z\^} is convergent, i.e. 
there is a unique zero limit measure. Here, by q{z) being symmetric we mean that after a translation on 
the real axis, q(z) is an even function. But for a non-symmetric quartic oscillator there are at least two 
zero limit measures. The Stokes lines play an important role in the description of the zero lines. In fact the 
infinite zero lines are asymptotic to Stokes lines. This fact was observed in _B . 

Our proofs are elementary. We use the complex WKB method, connection formulas and asymptotics of 
the eigenvalues by Fedoryuk in [FV . 

At the end of the paper (13. 2p we will briefly mention some interesting examples of one- well and double- well 
potentials where deg(g(z)) = 4, 6. 

Theorem 1.2. Let q{z) — (z^ - a^)(z^ - b"^), where < a < b. Then as n ^ oo 

Zx„ — > -\^/q{z)\ \d'-f\, 7 = (a, b) U (-6, -a) U {-ooi, +ooi). 

TT 

Notice in Theorem 11.21 we can express 7 by three equations 

7 = {US{a, z) = 0} U {nS{-a, z) 0} U {^S{a, z) = -i^}, 

where f = q{t) dt, and each equation is written in some canonical domain. 

Theorem 1.3. Let q{z) — (z^ - a^)(z^ + 6^), where a, 6 > 0. Then as n ^ 00 

Zx^ — » —\\/q(z)\ \d'^\, 7 = (—a, a) U {bi, +oo«) U {—ooi, —bi). 

TT 

14^6 also note that in Theorem \1.3\ we can express 7 by three equations 

7 = {^S{a, z) = 0} U {^S{bi, z) = 0} U {^S{-bi, z) = 0}, 
where each equation is written in some canonical domain. 

Theorems 11.21 and II .31 state that for a symmetric quartic polynomial there is a unique zero limit measure. 
This is not always the case when g(z) is not symmetric. Let q{z) — {z — ao)(z — ai)(z — a2){z — 03) where 
oo < ai < a2 < 03. Using the quantization formulas (for example in [S| IFlj ). we have two sequences of 
eigenvalues 

(7) AW = ^^,r + 0(l), ai = T \VW)\dt. 

2ai n Jao 

(8) A(2) = — ^7r + 0(-), a2= \VW)\dt. 

Now with this notation we have the following theorem: 

Theorem 1.4. Letq(z) = (z — ao)(z — ai)(z — a2)(z — 03), where ao < ai < 02 < 03 are real numbers. Then 

1. if ^ is irrational then for each £ G {lj2} there is a full density subsequence {Xnl} of {Xn^} such 
that 

^"k TT 

where 

71 = (flo, ai) U (02, as) U {5R(5(a2, z) = 0}, 

(9) 

72 = (ao, ai) U (a2, 03) U {^{S{ai,z) = 0}. 
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Figure 1. The zero lines of the non-symmetric quartic osciUator. The thick hnes are the 
zeros lines 71 and 72 given by ([9]). 

2. if — is rational and of the form , or '^^^^^ then for each £ € |1, 2| 

" TT 

3. if ^ is rational and of the form 2rl+i where gcd(2ri + 1, 2r2 + 1) ~ 1, then for each I G {1,2} there 
exists a subsequence {Xnl} of {Xn^} of density 2rl+i ■^^^/i that 

Z.w — -\^/q{z) \ 1^7^ I . 

-^"fc TT 

In fact {X'il} = {X[f^ \2n+l^Q {mod 2ri + 1)}. 

Figure ([T]) shows 71 and 72 defined in ([9]). As we see the zero lines here are some of the Stokes lines. We 
should mention that in Theorem II. 41 when — is irrational or of the form tP~^ 7^ 1, we do not know what 
happens to the rest of the subsequences. There might be some exceptional subsequences (of positive density 
in the case — — l^^^4) for which the zero lines are different from 7^. As we saw in Theorem II .21 this is the 
case for the symmetric double- well potential when ^ = 1. We probably need a more detailed analysis of 
the eigenvalues in order to answer this question. 

Remark 1.5. In this remark we mention some very recent results of A. Eremenko, A. Gabrielov, B. Shapiro 
in |EGS1| and |EGS2| . and compare them to ours as the interests and the approaches in these articles are 
very similar to ours. 

1. Theorems ] 1. Si and \1.3\ do not say anything about the exact location of the zeros but they only state that 
as n 00 the zeros approach to 7 with the distribution law in J^. It is easy to see that for both of these 
symmetric cases, for each n, all the zeros o/?/(z,A„) except finitely many of them are on 7. In |EGS1| . the 
authors prove that for the solutions of the equation 

(10) - y" + P{x)y = Xy, y E L^R), 

where P{x) is an even real monic polynomial of degree 4, the zeros of y{z) belong to the union of the 
real and imaginary axis. This result indeed implies that for all n, all the zeros of y{z, Xn) in Theorem 
and Theorem \1.3\ are on the corresponding 7. 

2. In }EGS2| . the authors show that the complex zeros of the scaled eigenfunctions F„(z) = y{xl/'^z) of 
l\l(M . where d =degree{P{x)), have a unique limit distribution in the complex plane as A„ — > cxd. The scaled 
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eigenf unctions satisfy an equation of the form 

r^'(z) = A:2(z'^-l + o(l))y„(z), 

The main reason that they could establish a uniqueness result for the limit distribution of complex zeros 
of Yn{z) is due the special structure of the Stokes graph of the polynomial z'^ — 1 which is proved in Theorem 
1 in |EGS2j . 

2. A Review of the Complex WKB Method 

To prove the theorems we first review some basic definitions and facts about complex WKB method. We 
follow [Fl| . See |0H IS] IEF] for more references on this subject. 



oo, 



We consider the equation 

(11) y"{z,\) = X'q{z)y(z,\), A 
on the complex plane C, where q{z) is a polynomial with simple zeros. 

2.1. Stokes lines and Stokes graphs. A zero zqj of q{z) is called a turning point. We let S{zq,z) = 
J^^ \J q{t)dt. This function, is in general, a multi-valued function. The maximal connected component of the 
level curve ^{S{zo, z)) — with initial point zq and having no other turning points are called the Stokes 
lines starting from zq. Stokes hnes are independent of the choice of the branches for S{zo, z). The union of 
the Stokes lines of all the turning points is called the Stokes graph of (fTTj) . 

Figure ([2]) shows the Stokes graphs of many polynomials. 

Since the turning points are simple, from each turning point three Stokes lines emanate with equal angles. 
In general if zo is a turning point of order n, then n + 2 Stokes lines with equal angles emanate from zq- 

2.2. Canonical Domains, Asymptotic Expansions. Since q{z) is a polynomial, the Stokes graph divides 
the complex plane into two type of domains: 

1. Half-plane type: A simply connected domain D which is bounded by Stokes lines is a half-plane type 
domain if under the map S = S{zo, z), it is biholomorphic to a half-plane of the form KS* > a or 

< a. Here zq is a turning point on the boundary of D. 

2. Band-type: D as above, is of band-type if under S, it is biholomorphic to a band of the form 
a<nS <b. 

A domain D in the complex plane is called canonical if S{zo, z) is a one-to-one map of D onto the whole 
complex plane with finitely many vertical cuts such that none cross the real axis. A canonical domain is the 
union of two half-plane type domains and some band type domains. For example, the union of two half-plane 
type domains sharing a Stokes line is a canonical domain. 

Let e > be arbitrary. We denote D'^ for the pre-image of S{D) with e-neighborhoods of the cuts and 
£-neighborhoods of the turning points removed. A canonical path in D is a path such that K(S') is monotone 
along the path. For example, the anti-Stokes lines (lines where, 3(5*) = 0), are canonical paths. For every 
point z in , there are always canonical paths J~{z) and 7^(2;) from z to 00, such that SftS J, —00 and 
5R5 t 00, respectively. 
Now we have the following fact: 

With D, 7+, 7~ as above, to within a multiple of a constant, equation ([TT|) has a unique solution yi{z, A), 
such that 

(12) lim yi(z,A) = 0, SR5(zo, z) i -00, 

and a unique solution i/2(z, A) (up to a constant multiple), such that 

(13) lim 2/2(^1 A) =0, RS{zo, z) 00. 

The solutions yi and 2/2 in (|12p and (|13p have uniform asymptotic expansions in in powers of j. Here 
we only state the principle terms: 



(14) yi(z,A)-gT^(z)e^^(^«^^)(l + ei(z,A)) A 



00, 
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q(z) = (z^-a^)(z^+bz+c) q(z) = (z-a(,)(z-a ,)(z-a,)(z-a3 ) 



Figure 2. Stokes lines for some polynomials 

(15) 2/2(z,A) = gT(^)e-^5(.o-)(l^_£2(z,A)) A ^ oo, 
where 

(16) ei{z, X) — 0{ — ), uniformly in Z)"^ , A — > oo, 

A 

(17) £2(2, A) = 0( y), uniformly in Z?"^ , A ^ 00. 

A 

Notice that the equalities (|16p and ([T7]) would not necessarily be uniformly in D' if q(z) was not a 
polynomial. 

2.3. Elementary Basis. Let D be a canonical domain, / a Stokes line in D, and zq G I a turning point. We 
use the triple {D,1,zq) to denote this data. We select the branch of S{zo,z) in D such that ^S{zo,z) > 
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for z I. The elementary basis {u{z),v{z)} associated to {D,1,zq) is uniquely defined by 

A) = cyi(z, A), and v{z, X) ^ cy2{z, X), 

\c\ = 1, and arg(c) = lim^^^^.^ei a'rg{q'^/*{z)), 
where yi{z,X), y2{z,X) are given by IHl) and 



(18) 



2.4. Transition Matrices. Assume {D,l,zo)j and {D,l,zo)k are two triples and = {uj,Vj} and = 
{wfejWfc} their corresponding elementary basis. The matrix Qjk{X) which changes the basis Pj to f3k is called 
the transition matrix from (3j to (3k- 

Fedoryuk, in [EF], introduced three types of transition matrices that he called elementary transition 
matrices, and he proved that any transition matrix is a product of a finitely many of these elementary 
matrices. The three types are 

1) {D,l, zi) i—f (D,l, Z2). This is the transition from one turning point to another along a finite Stokes 
line remaining in the same canonical domain D. The transition matrix is given by 

(19) n{X} = e^^ ^1 '^l'" y \S{z,,Z2)\, e'^ = ^. 

2) {D, li, zi) ^ {D, I2, 22)- Here the rays S{li) and S{l2) are directed to one side. This is the transition 
from one turning point to another along an anti-Stokes line, remaining in the same domain D. The 
transition matrix is 

(20) n{X) = e'^ ""p'" ) , a = \S{z,,Z2)\, e'^ = |. 

3) (I?i,Zi,zo) ^ {D2,l2, zq) This is a simple rotation around a turning point zq so that Di and D2 
have a common sub-domain. More precisely, let = 1,2,3} be the Stokes lines starting at zq 
and ordered counter-clockwise so that Ij+i is located on the left side of Ij. We choose the canonical 
domain Dj so that the part of Dj on the left of Ij equals the part of D^+i on the right of Ij+i- Then 



(21) 



"^'^ I 1 ^a,+i.+2(A) j ' 

a,,,+i(A) = l + 0(i), l<j<3, 

ai,2(A)Q!2,3(A)Q!3,i(^) = 1; and ajj+i(A)aj+ij (A) = 1. 



2.5. Polynomials with real coefficients. We finish this section with a review of some properties of the 
Stokes lines and transition matrices in ([2T|) when the polynomial q{z) has real coefficients. 

1) The turning points and Stokes lines are symmetric about the real axis. If xi < X2 are two real 
turning points and q{x) < on the fine segment / = [xl,X2]^ then / is a Stokes fine (See Figure 
Similarly, if q{x) > on Z, then / is an anti-Stokes line. Let xq be a simple turning point on the real 
axis, and let Iq, 11,12 be the Stokes lines starting at xq. Then one of the Stokes lines, say Iq, is an 
interval of the real axis, and I2 — h- The Stokes lines h and I2 do not intersect the real axis other 
than at the point xo- If a Stokes line I intersects the real axis at a non-turning point, then I is a 
finite Stokes line and it is symmetric about the real axis. 

If Imix^oo q{x) — 00, and a;+ is the the largest zero of q{x), and lo,h,l2 are the corresponding 
Stokes lines, then there is a half type domain Z)^ such that 

[x+,+oo]cD+, D+=D+, liUl2CdD+. 

Clearly [a:o,+oo] is an anti-Stokes line and S{xq,oo) = 00. By (|12p . there exists a unique solution 
y+(z. A) such that 

lim y^{x, A) = 0. 
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(22) 



(23) 



Figure 3. 

Similarly by (|13p if lim3.^_oo 9(2;) = 00 and x~ is the smallest root of q{x), and D~ a half type 
domain containing [— oo,xo], there exists a unique solution y~(z,X) such that 

lim y^{x, A) = 0. 



X — ^ — 00 



Therefore if y{x, A) is an L^-solution to ([2]) , then for some constants 

y{x, A) = c+y+{x, A) = c^y^{x, A). 
Now let ri_|-__(A) be the transition matrix connecting D+ to D'^ and let 



The fact that y'^{x, A) is a constant multiple oi y (x, A) is equivalent to 

b{\) = 0, 

which is the equation that determines the eigenvalues A„. To calculate 1]+__(A) and hence 6(A) 
we have to write this matrix as a product of finitely many elementary transition matrices connecting 
D+ to D-. 

2) When the polynomial q{z) has real coefficients, the transitions matrices in (|21l) have some symmetries. 
Let Xq be a simple turning point and q{x) > on the interval [xQ,b]. We index the Stokes lines 
lo,li, h as in Figure We define the canonical domains Dq, Di, D2 by their internal Stokes lines 
and their boundary Stokes lines as the following 

Do^D^, loCDoJiUhcdD, 

[xo,b]cDi, loUhcdDi, 

Now with the same notation as in (PT|) . we have 

1. 



"04 = "0,2, |ai,2 

3. Proofs of the Theorems 

First of all we have the following lemma: 

Lemma 3.1. Let T = (D, I, zq) be a triple as in section 2.3 and let {u{z, A), v{z, A)} be the elementary basis 
associated to T in hll^). We write y{z, A„) in this basis as 

(24) y{z, A„) = a(A„)u(z, A„) + b{Xn)v{z, A„). 

//{A„^. } is a subsequence o/{A„} such that the limit 



(25) 



t 



lim 



1 



k — >co 2A7J 



■log 



a(A„J ' 
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exists, then in we have 

^A„, ^-Iv^IMtI, -f^{zeD\^Sizo,z)^t}. 

The last expression means that for every ip e C^{D^) we have 



Zx^S'fi) ^ - / V{z)\^q(z)\\d-f\. 

Proof of Lemma. For simplicity we omit the subscript in A„j. , but we remember that the hmit in 
(^5)) is taken along X^^. Using ([M]) . and the equation y{z,X) = in 13^, is equivalent to 

(26) 5(zo, z) - - log( Y^^) = 2^ log I ^ I +.(^. + 2X -g(^)), e 

where we have chosen log z ~ \ogr + i9, — tt < 9 < tt. We use S{z) for the function on the left hand side 
of (|26p and for the sequence of complex numbers on the right hand side. As we see S{z) is the sum of the 
biholomorphic function S{zq, z) and the function 

(27) ii{z, A) -i- log(i±i4^) = O(^), uniformly in by (HZD. 

2A l + £2(-z, A) A^ 

Now suppose ¥5 e Cf[D^) and = supp((^). We also define K' = S{K) where S{z) = S{zo, z). Without 
loss of generality we can assume that {x — t}r\int{K') is a connected subset of the vertical line x — t, because 
we can follow the same argument for each connected component. Now let s — lengthdx — t} f) int(iir')). It 
is clear that because of (1251) 



(28) N := #{ flfc G int(X')} sA. 

We call this finite set {ak}m+i<k<m+N ■ Now let X C 1^ C Z?^ be an open set with compact closure in 
D^. We choose A large enough such that 

|Ai(z,A)| < |S'(zo,z) -a| yaeK', y z e dV. 

Since is a biholomorphic map, by Rouche's theorem the equation 

S{z) = ttk, m + l<k<m + N, 

has a unique solution Zk in V for each k. Now by (P5)) . (pS)) . and (P7)) . we have 



9k + ^ 1 
(t + 0(1) + ^(^^^ - 3(M(^fe, A)) + 0(-))) . 



It follows that 



m+N , m+N 91 M ^ 1 

^a(¥') = ^ E = X ^ (^°5-i)(i + o(l) + *(^^^-S(M(^fe,A)) + 0(-))). 

fc— m+l k—m-\-l 

Using the mean value theorem on the a;— axis and ([28]), we obtain 
(29) lim Z,((^) = hm - ^ {[(^ o S-'){t + ^{—^n - ^{^^{z,, X)) + 0{-)))] + o{l)}. 

A — *oo A — >OQ A ^"^^ LA A 

Because of ([27|l . we know that 0^(^(2;^, A)) = 0{^) uniformly in k. Therefore the set 

P = {(t, ^^TT - S(M(^fe, A)) + 0(-))| m + 1 < < m + TV} 

is a partition of the vertical interval \x = nint(i^') with mesh(p) ^ as A ^ oo. This together with 
implies that 
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lim Zx{(p) — — If o S ^dy. 

TT J{x=t} 

Now, if in the last integral we apply the change of variable z ^ S{z), then by the Cauchy-Riemann 
equations for S, we obtain 

-/ ipoS^^dy^- ( ip{z)\^/l^\\d-i\. 

This proves the Lemma. 

Proof of Theorem 11.11 First of all, we cover the plane by finitely many canonical domains Dm- Let 
£ > be sufficiently small as before. Assume {Z\^^ } is a weak* convergent subsequence converging to a 
measure Z. Clearly {Z\^^ } converges to Z in each Df^^. We claim that the limit exists for every triple 
Tm = {Dm, Zm,lm)- This is clcar from Lemma 13.11 This is because if in (j25p we get two distinct limits 
ti and t2 for two subsequences of {Z\^^}, then we get two corresponding distinct limits Zi and Z2 which 
contradicts our assumption about {Z\^^ }. We should also notice that if in (pS)) . t ~ +00 then in the proof 

of Lemma [01 for A large enough we have {x = log | |} fl mt{K') = 0, and therefore Z\]je^ = 0. This 
means that we do not obtain any zero lines in this canonical domain. In other words the zeros run away from 
this canonical domain as A — > 00. But as we mentioned in the introduction, all the Stokes lines on the real 
axis are contained in the set of zero lines of every limit Z, meaning that in Theorem ll.il 7 is never empty. 

Now notice that because IJ^ Dm covers the plane except the e-neighborhoods around the turning points, 
we have proved that 

Z{^) - i /, V>{Z)\V^)\ M7I, ^ G C-(C\ Dm B{Zrn, s)) . 

To finish the proof we have to show that if ip^ G C^{[J^ B{zm,£)) is a bounded function of e, then 

lim lim sup (ipe) — 0. 
This is clearly equivalent to showing that if zq is a turning point, then 



(30) lim limsup #{^ ^ ^(^0. f A) ^ 0} ^ ^ 

To prove this we use the following fact in [Fl] pages 104 — 105 or [EFj pages 39 — 41, which enables us 
to improve the domain Z?^ in p6)) and (fT7|) from a fixed e to e{X) dependent of A such that e(A) ^ as 
X CO. 

Let -D be a canonical domain with turning points Zm on its boundary. Assume N{X) is a positive function 
such that N{oo) = 00. Now if we denote 

DiX) = D \\jB{zm, k'(z„)ri/37V(A)A-2/3), 

rn 

Then in place of equations ((T6)) and ((T7|) we have 

ei{z,X), e2iz,X) = 0{N{X)~^/'^), uniformly in !?( A), A ^ 00. 

In fact this implies that Lemma l3.1l is true for every ip supported in D. This is because we can follow the 
proof of the lemma line by line except that in (^7)1 we get A) — 0{N{X)^^^^ X^^) uniformly in D{X) and 
therefore, using N{oo) = 00, we can still conclude mesh(p) — > as A — > 00. 

We choose A^(A) = X^^^'^. By the discussion in the last paragraph in ([50)) we can replace e by e(A) = 
cN{X)X^^^'^ = cA^^/^^, where c = |'?'(-Zo)l"^^'^- Let us find a bound for the number of zeros of y{z,X) in 
B{zo,e{Xj). Let M = snp^j^^^ g^{\q{z)\) where 5 > is fixed and is chosen such that the ball B{zo,d) does 
not contain any other turning points. We also choose A large enough so that e(A) < 5. If C is a zero of y{z, A) 
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Figure 4. 



in the ball B{zo,e{X)) then by Corollary 11.1.1 page 579 of |H1| we know that there are no zeros of y{z, A) 

herefore 

area(B(zo,e(A) + ^A-i)) 



in the ball of radius -^^A ^ around ( except (. Therefore 



#{^6i?(zo,£(A))|y(z,A) = 0}< 



area(i?(C,^A-i)) 



= 0(A'^/«), 



and so 



lim 

A — >oo 



#{zei3(zo,e(A))|y(z,A)-0} 



This finishes the proof of Theorem 11.11 



A 



The Proofs of Theorems 11.21 11.41 We wiU not prove Theorem II. 3( because the proof is similar to 
(in fact easier than) the proof of the two well potential. To simplify our notations let us rename the turning 
points as xi = oq, Xm = o,i,Xn — 0-2, Xp = 03. Then we can index the Stokes lines as in Fig ([?]). 

We define the canonical domains Di, Dmo, D, Dno, and Dp by 



h C Di, mi,mo,l2 C dDi, 

mo C Dmo, hj2,mi,m2 C dDmo, 
(31) mi,niCD, li,lo,m2,n2,PQ,Pi C dD, 

noCDno, Pi,P2,ni,n2 C dDno, 
Pi C Dp, ni,no,p2 C dDp . 

Notice that the complex conjugates of the these canonical domains are also canonical domains and in fact 
if we include these complex conjugates then we obtain a covering of the plane by canonical domains. But 
because q(x) is real, the zeros are symmetric with respect to the a;-axis, and it is therefore enough to find 
the zeros in Di U Dmo U D U U Dp. By lemma [3T] we only need to discuss the limit (|25l) in each of these 
canonical domains. First of all let us compute the equation of the eigenvalues ([22]) . 

Here, the transition matrix is the product of the seven elementary matrices associated to the fol- 

lowing sequence of triples: 

{Dp,li,Xp) 1-^ {Dno,nQ,Xp) 1-^ (D„(,,no,a;„) 1-^ {D,ni,x„) 

^{D ) ^ {Dmo,mo,Xm) {Dmo,mo,Xl) (Di,li,Xl). 
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In fact if we define 

ai=/ W9it)\dt, «2=/ W9it}\dt, / V9it)dt, 

J xi 

then by (ITn]),(12DD and (HI]), we have 



6(A) 1 



«;:on. ^ r e-^"^ \ / a-/p„ A / 



A in I \ p*-^"2 n / V 1 70' / V 1 

A simple calculation shows that 

Hence 6(A) = implies that 

(32) ri(A)r2(A) = a,/,„a,;„i„^e'^(— ^)e-2^«, 

where 

ri(A) = a„,,™,e-*^"^ + a,,i,a„i„^e^^"i = 2cos(aiA) + 0(i), 
^^^^ r2(A) = ap„p,e-*^"^ + a„,„,a-ip„e*^"^ = 2cos(a2A) + 0[\). 

Now let us discuss the limit in ((25|) for each of the canonical domains defined in (|3T|) . Even though the 
coefficients ci(A), 6(A) are different for different canonical domains, wc do not consider it in our notation. 
By (fT4|) . (fT5|) , (fT6|) . and (fT7|) . it is clear that for A large enough there are no zeros in Df and D^. For 

(L'„o,no,a:p) we have 

«(A) A / a-ip„ U M = <vo 
6(A) y I 1 zapop2 y V I 



Using ipT|) . ([^ . for the full sequence A„ we have 

1 1 |&(An)| 1,1 I 

2A;^°si^i-23ri°si"--i 



Hence t = and, by Lemma l3.H the Stokes line no = {0,2, 03) is a zero line in Duq- The same proof shows 
that the Stokes line [ao,ai] is a zero line for the full sequence in Dmo- Now it only remains to discuss 
the limit in (pS)) in the canonical domain D. For the triple (I?,ni,x„) we have 

a(A) \ / «-„i„^ \ / e-»^- \ / a-^^^ \( = ( ^''"'<Po<n, 
6(A) ; 1^ 1 m„,„, ; 1^ e^^"= J\l iap„p, J\l J \ ^T2{\) 



Therefore, by the second equation in (I2I|) . we obtain 

(34) lim ^log| I = lim log |r2(A„)|. 

The limit ([M)) does not necessarily exist for the full sequence {A„}. We study this limit in different cases 
as follows: 
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This is exactly the symmetric case in Theorem ll.2l It is easy to see that if ai = a2, then there exists 
a translation on the real line which changes q{z) to an even function. When q{z) is even, because of 
the symmetry in the problem, we have ri(A) = r2(A). On the other hand equation (j32p implies that 

|ri(A)||r2(A)He-2^«(i + o(l)). 

This means that in the symmetric case, the full sequence A„ satisfies 

|ri(A„)| = |r2(A„)|=e-^"«(l + 0(^)). 

Therefore, by (I34p wc have i = — and using the lemma the line diS{a, z) — — is the zero line 
in D. We note that this in fact determines the whole imaginary axis, because diS{a, 0) — — ^C- Also 
notice that in the symmetric case, by our notations we have a — a2 = Xn. This proves Theorem ll.21 2. 

Ot-2 ' 

In this case as we mentioned in the introduction, there are more than one zero limit measures. Here 
the limit ([M]) behaves differently for the two subsequences in ([7]) and ([5]) (notice that the equations 
([7]) and ([8]) in fact follow from ([32])). It is clear from ((33|) that if for a subsequence {A„^, } we 
have a lower bound 8 for | cos(q:2A„^)|, then we have t^limfe^oo 2X~ l-'^2(A„j. )| — 0. Also if we 

have a lower bound 8 for |cos(aiA„^)| then limfe^oo 2X~~ ^'-'S l-ri(A„^. )| = 0, and by (|35p we have 

t ~ linifc^oo 2\ ^"-"S |r2(A„^)| — To find such subsequences we denote for each £ = 1,2 

Af^ = {A„; |cos(a£A„)| > 5}. 

By O and ([5]), it is clear that up to some finite sets A^^^ C {Xn^} and Af^ C {A^^'}. We would 
like to find the density of the subsets A'^'^ and Af^ in {Xn ^} and {A^^-*} respectively. Here by the 
density of a subsequence {A„^} of {A„} we mean 

hm^i^iAi^^. 



If we set r = arcsin((5) then we have 

h- + irr^ + b\>r + 0{- 
2 a2 2 n 



= e N; \{n + -)— + {m+^)\>T + 0{^), V m £ Z}, 



^ £ |(„ + i)^ + + i)| > ^ + o{-), V m e Z}. 
2 Qfi 2 n 

We only discuss the density of the subset A^'^ . We rewrite this subset as 

= {n e N; \(2n + l)ai + (2m + 1)02! > 2a2T + O(-), V m G Z}. 

n 

From this we see that if ^ is a rational of the from 2r2+i ^'^^ then because for every m and 

n we have 

\{2n + l)(2ri) + (2m + l)(2r2 + 1))| > 1, 
therefore = 1 for t = This proves Theorem ll.41 2. When — is a rational of the form 

we define 

S^^^ ={ne 2n+1^0 (mod 2r2 + 1)}. 

Since for every n G B^g'^ and m G Z we have 

\(2n + l)(2ri + 1) + (2m + l)(2r2 + 1))| > 1, 
for T = g;:ip4, we get d(4^^) > '^(^f ^) = 2^^- This completes the proof of Theorem [UlS. 



14 HAMID HEZARI 




q{z) = (z2 - a2)(z2 + bz + c) q{z) = {z^ - a^){z^ - b'')[z'' + c^) 



Figure 5. 

To prove Theorem 11.41 1. when ^ is irrational, we use the fact that the set Zai © Za2 is dense in 
R. In fact it is easy to see that the subset A = {nai + raa2\ n £ N, m e Z} is also dense. Now if 
we rewrite A^^ as 

Af'' = {new, \{nai + maa) + ^(ai + a2)| > + O(^), V m e Z}, 

then from the denseness of the set A, it is not hard to see that in this case d(A^p) = 1 — ^^r. 

Hence we conclude that when / = 1 there is a subsequence {An^^} of {An''} of density 1. The same 
argument works for 1 = 2. This finishes the proof. 

Remark 3.2. In Figure ^ we have illustrated the zero lines for the polynomials 

q{z) = {z'^ ^a^){z^ + bz + c) and q{z) = {z"^ - a^){z^ ~ b^){z^ + c^) . 

The thickest lines in these figures are the zero lines. In fact for these examples there is a unique zero limit 
measure as in the other symmetric cases we mentioned in Theorems and \l . We will not give the proofs, 
as they follow similarly, but we would like to ask the following question: 

Is there any polynomial potential with n wells, n > 3 for which there is a unique zero limit measure for 
the zeros of eigenf unctions ? 

Acknowledgements: I am sincerely grateful to Steve Zelditch for introducing the problem and many 
helpful discussions and suggestions on the subject. 
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